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We present an extensive analysis of squark and gaugino hadroproduction and decays in non- 
minimal flavour violating supersymmetry. We employ the so-called super-CKM basis to define the 
possible misalignment of quark and squark rotations, and we use generalized (possibly complex) 
charges to define the mutual couplings of (s)quarks and gauge bosons/gauginos. The cross sections 
for all squark-(anti-)squark/gaugino pair and squark-gaugino associated production processes as 
well as their decay widths are then given in compact analytic form. For four different constrained 
supersymmetry breaking models with non-minimal flavour violation in the second/third generation 
squark sector only, we establish the parameter space regions allowed/favoured by low-energy, elec- 
troweak precision, and cosmological constraints and display the chirality and flavour decomposition 
of all up- and down-type squark mass eigenstates. Finally, we compute numerically the dependence 
of a representative sample of production cross sections at the LHC on the off-diagonal mass matrix 
elements in the experimentally allowed/favoured ranges. 
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INTRODUCTION 



Weak scale supersymmetry (SUSY) remains a both theoretically and phenomenologically attractive extension of 
the Standard Model (SM) of particle physics Apart from linking bosons with fermions and unifying internal and 

external (space-time) symmetries, SUSY allows for a stabilization of the gap between the Planck and the electroweak 
scale and for gauge coupling unification at high energies. It appears naturally in string theories, includes gravity, and 
contains a stable lightest SUSY particle (LSP) as a dark matter candidate. Spin partners of the SM particles have not 
qq , yet been observed, and in order to remain a viable solution to the hierarchy problem, SUSY must be broken at low 
' energy via soft mass terms in the Lagrangian. As a consequence, the SUSY particles must be massive in comparison 
to their SM counterparts, and the Tevatron and the LHC will perform a conclusive search covering a wide range of 
f~>) i masses up to the TeV scale. 

[~^. | If SUSY particles exist, they should also appear in virtual particle loops and affect low-energy and electroweak 
precision observables. In particular, flavour-changing neutral currents (FCNC), which appear only at the one-loop 
level even in the SM, put severe constraints on new physics contributions appearing at the same perturbative order. 
• i-h , Extended technicolour and many composite models have thus been ruled out, while the Minimal Supersymmetric 
^ ■ Standard Model (MSSM) has passed these crucial tests. This is largely due to the assumption of constrained Minimal 
Flavour Violation (cMFV) [1, [| or Minimal Flavour Violation (MFV) @, H, 0, where heavy SUSY particles may 
appear in the loops, but flavour changes are either neglected or completely dictated by the structure of the Yukawa 
couplings and thus the CKM-matrix [I, Q . 

The squark mass matrices M?, Mj~, and M|> are usually expressed in the super-CKM flavour basis In MFV 

SUSY scenarios, their flavour violating non-diagonal entries Ay, where i,j — L,R refer to the helicity of the (SM 
partner of the) squark, stem from the trilinear Yukawa couplings of the fermion and Higgs supermultiplets and the 
resulting different renormalizations of the quark and squark mass matrices, which induce additional flavour violation 
at the weak scale through renormalization group running [lj [l2|, [H, EH, while in cMFV scenarios, these flavour 
violating off-diagonal entries are simply neglected at both the SUSY-breaking and the weak scale. 

When SUSY is embedded in larger structures such as grand unified theories (GUTs), new sources of flavour violation 
can appear [l5| . For example, local gauge symmetry allows for i?-parity violating terms in the SUSY Lagrangian, but 
these terms are today severely constrained by proton decay and collider searches. In non-minimal flavour violating 
(NMFV) SUSY, additional sources of flavour violation are included in the mass matrices at the weak scale, and 
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their flavour violating off-diagonal terms cannot be simply deduced from the CKM matrix alone. NMFV is then 
conveniently parameterized in the super-CKM basis by considering them as free parameters. The scaling of these 
entries with the SUSY-breaking scale Msusy implies a hierarchy All 3> Alr.rl ^> Arr [La ]. 

Squark mixing is expected to be largest for the second and third generations due to the large Yukawa couplings 
involved [16j |. In addition, stringent experimental constraints for the first generation are imposed by precise mea- 
surements of K° — K° mixing and first evidence of D° — D° mixing [I?], E^l • Furthermore, direct searches of 
flavour violation depend on the possibility of flavour tagging, established experimentally only for heavy flavours. We 
therefore consider here only flavour mixings of second- and third-generation squarks. 

The direct search for SUSY particles constitutes a major physics goal of present (Tevatron) and future (LHC) 
hadron colliders. SUSY particle hadroproduction and decay has therefore been studied in detail theoretically. Next- 
to-leading order (NLO) SUSY-QCD calculations exist for the production of squarks a nd g luinos (20j . sleptons [2l| . 
and gauginos [12] as well as for their associated production [23| . The production of top [24| and bottom [25| squarks 
with large helicity mixing has received particular attention. Recently, both QCD one-loop and electroweak tree-level 
contributions have been calculated for non-diagonal, diagonal and mixed top and bottom squark pair production [2(| . 
However, flavour violation has never been considered in the context of collider searches apart from the CKM-matrix 
appearing in the electroweak stop-sbottom production channel [261 ] . 

It is the aim of this paper to investigate for the first time the possible effects of non-minimal flavour violation 
at hadron colliders. To this end, we re-calculate all squark and gaugino production and decay helicity amplitudes, 
keeping at the same time the CKM-matrix and the quark masses to account for non-diagonal charged-current gaugino 
and Higgsino Yukawa interactions, and generalizing the two-dimensional helicity mixing matrices, often assumed to be 
real, to generally complex six-dimensional helicity and generational mixing matrices. We keep the notation compact 
by presenting all analytical expressions in terms of generalized couplings. In order to obtain numerical predictions 
for hadron colliders, we have implemented all our results in a flexible computer program. In our phenomenological 
analysis of NMFV squark and gaugino production, we concentrate on the LHC due to its larger centre-of-mass energy 
and luminosity. We pay particular attention to the interesting interplay of parton density functions (PDFs), which 
are dominated by light quarks, strong gluino contributions, which are generally larger than electroweak contributions 
and need not be flavour-diagonal, and the appearance of third-generation squarks in the final state, which arc easily 
identified experimentally and generally lighter than first- and second-generation squarks. 

After reviewing the MSSM with NMFV and setting up our notation in Sec. [H] we define in Sec. IIIII generalized 
couplings of quarks, squarks, gauge bosons, and gauginos. We then use these couplings to present our analytical 
calculations in concise form. In particular, we have computed partonic cross sections for NMFV squark-antisquark 
and squark-squark pair production, squark and gaugino associated and gaugino pair production as well as NMFV 
two-body decay widths of all squarks and gauginos. Section IIVI is devoted to a precise numerical analysis of the 
experimentally allowed NMFV SUSY parameter space, an investigation of the corresponding helicity and flavour 
decomposition of the up- and down-type squarks, and the definition of four collider-friendly benchmark points. These 
points are then investigated in detail in Sec. [V] so as to determine the possible sensitivity of the LHC experiments on 
the allowed NMFV parameter regions in the above-mentioned production channels. Our conclusions are presented in 
SecLYIl 



II. NON-MINIMAL FLAVOUR VIOLATION IN THE MSSM 



Within the SM, the only source of flavour violation arises through the rotation of the up-type (down-type) quark 
interaction eigenstates u' L R (d' L R ) to the basis of physical mass eigenstates ul,r (rfi.fl), such that 

d L ,R = V£ R d' LR and u l ,r = V[ l R u' L R . (1) 

The four bi-unitary matrices V^' R diagonalize the quark Yukawa matrices and render the charged-current interactions 
proportional to the unitary CKM-matrix 0, Q 

(V u d v us V u b \ 
Vcd V cs Vcb ■ (2) 
Vd V ts V tb J 

In the super-CKM basis [lOj], the squark interaction eigenstates undergo the same rotations at high energy scale 
as their quark counterparts, so that their charged-current interactions are also proportional to the SM CKM-matrix. 
However, different renormalizations of quarks and squarks introduce a mismatch of quark and squark field rotations 
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at low energies, so that the squark mass matrices effectively become non-diagonal [TT|, d, EH 13 ■ NMFV is then 
conveniently parameterized by non-diagonal entries Af? with i,j = L, R in the squared squark mass matrices 



M- = 



2 

A?" 



LL 



A MC 

M l 

A tc C 



and 



m u X u 

A CM 

LR. 

A til 

^LR 



( M\ 
All 
A£1 



A tic 

^RL 

m c X c 

Ate 
^LR 



Ads 

M l 
Abs 3 



Aut 
Act 

M? 



\ 



m d X d 

Asd 
^LR. 

Abd 
^LR 



Ads 
^RL 

m s X s 

Abs 
^LR 



^RL 

Act 

^RL 

m t X t 



Adb 
Asb 

M~ 



m u X u 

ACU 

^RL 

Atu 

^RL 



A tic 

^LR 

m c X c 

Ate 

^RL 



m t X t 



&RR 



AUC 

^RR 

M l 

Ate 

^RR 



^■RR 

Act 

^RR 



(3) 



A* 
^RL 

Asb 
^RL 

m b X b 



m d X d 


Ads 
^LR, 


Adb 
^LR 




Asd 
^RL 


m s X s 


Asb 
^LR 




Abd 
^RL 


Abs 
^RL 


m b X b 




M k 


Ads 
^RR 


Adb 
^RR 




Asd 
^RR 


M l 


Asb 
^RR 




Abd 
^RR 


Abs 
^RR 


M k 


) 



(4) 



Mf = M% 



where the diagonal elements are given by 

-cos2/3M|(T g 3 -e 9 s^), 
cos 2/3 M^e q s^Y for up — type squarks, 
;/ ,- cos2/3M§e g s^ for down — type squarks, 
while the well-known squark helicity mixing is generated by the elements 
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cot (3 for up — type squarks, 
tan (3 for down — type squarks. 



(5) 
(6) 
(7) 

(8) 



Here, m q , T q , and e q denote the mass, weak isospin quantum number, and electric charge of the quark q. raz is 
the Z-boson mass, and sw [cw) is the sine (cosine) of the electroweak mixing angle 8w The soft SUSY-breaking 
mass terms are Mq and Mjy ^ for the left- and right-handed squarks. A q and /i are the trilinear coupling and 
off-diagonal Higgs mass parameter, respectively, and tan/3 = v u /v d is the ratio of vacuum expectation values of the 
two Higgs doublets. The scaling of the flavour violating entries AJJ with the SUSY-breaking scale Msusy implies a 
hierarchy A^ 3> Al R rl 3> A RR among them [l5[ . They are usually normalized to the diagonal entries [lH , so that 



A" 



(9) 



Note also that SU{2) gauge invariance relates the (numerically largest) A qq L of up- and down-type quarks through 
the CKM-matrix, implying that a large difference between them is not allowed. 

The diagonalization of the mass matrices M? and Mj requires the introduction of two additional 6x6 matrices 
R u and R d with 



diag (m^ , 
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diag (ml,..., ml) = R d MjR d l (10) 
; m qe . The physical mass eigenstates are given by 



(11) 
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In the limit of vanishing off-diagonal parameters A** , the matrices R q become flavour-diagonal, leaving only the 
well-known helicity mixing already present in cMFV. 



4 



III. ANALYTICAL RESULTS FOR PRODUCTION CROSS SECTIONS AND DECAY WIDTHS 



In this section, we introduce concise definitions of generalized strong and electroweak couplings in NMFV SUSY 
and compute analytically the corresponding partonic cross sections for squark and gaugino production as well as their 
decay widths. The cross sections of the production processes 



ah a (Pa)b hb {Pb) 



(12) 



are presented for definite hclicitics h a ^ of the initial partons a, b = q,q, g and expressed in terms of the squark, 
chargino, neutralino, and gluino masses rriq k , w~±, m^o, and m g , the conventional Mandelstam variables, 

s = (p a +Pb) 2 , t = {p a - pi) 2 , and u = (p a - p 2 ) 2 ', (13) 

and the masses of the neutral and charged electroweak gauge bosons mz and mw- Propagators appear as mass- 
subtracted Mandelstam variables, 
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(14) 



x <3i - •■~g i 

Unpolarized cross sections, averaged over initial spins, can easily be derived from the expression 

da = ^ — , (15) 

while single- and double-polarized cross sections, including the same average factor for initial spins, are given by 
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A. Generalized Strong and Electroweak Couplings in NMFV SUSY 



Considering the strong interaction first, it is well known that the interaction of quarks, squarks, and gluinos, whose 
coupling is normally just given by g s = \/Aira s , can in general lead to flavour violation in the left- and right-handed 
sectors through non-diagonal entries in the matrices R q , 

{ L q '/• '/• R q ■! '/ ! = { R< jk'- R j(k+3)} ■ ( 18 ) 

Of course, the involved quark and squark both have to be up- or down-type, since the gluino is electrically neutral. 

For the electroweak interaction, we define the square of the weak coupling constant gyy = e 2 / sin 2 9w m terms 
of the electromagnetic fine structure constant a = e 2 /(47r) and the squared sine of the electroweak mixing angle 
xw = sin 2 9w — Syy = 1 — cos 2 9w — 1 — . Following the standard notation, the W ± — xl — xf > Z ~ xt — Xj > an d 
Z — Xi — Xj interaction vertices are proportional to 
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In NMFV, the coupling strengths of left- and right-handed (s)quarks to the electroweak gauge bosons are given by 

{Lqq'Z, Rqq'z} = (2 T q — 2 e q Xw) X Sqq' , 

3 

{Lq^Z^Rqiiiiz} = - 2eqX w ) X y^{fi" fc Rjfl, Rj( 3+k ) Rj*3 + k)}' 
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{Lqq'W, Rqq'w} = {V^Cff V,,',0}, 



(20) 



k,l=l 



where the weak isospin quantum numbers are T? = ±1/2 for left-handed and T 3 ? ~j = for right-handed up- and 
down-type (s)quarks, their fractional electromagnetic charges are denoted by e/gg\, and Vjy are the elements of the 
CKM-matrix defined in Eq. (|2|). To simplify the notation, we have introduced flavour indices in the latter, d\ = d, 
d 2 = s, g?3 = b, u± = u, u 2 = c, and 113 = t. 

The SUSY counterparts of these vertices correspond to the quark-squark-gaugino couplings, 
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(21) 



where the matrices N, U and V relate to the gaugino/Higgsino mixing (see App. All other couplings vanish due 
to (electromagnetic) charge conservation (e.g. L-. ui -±). These general expressions can be simplified by neglecting the 
Yukawa couplings except for the one of the top quark, whose mass is not small compared to mw- For the sake of 
simplicity, we use the generic notation 



{^abci^abc} ~ {Labc, Rabc} 



(22) 



in the following. 



B. Squark-Antisquark Pair Production 



The production of charged squark-antisquark pairs 

q(h a ,p a )cf(h b ,p b ) -> Ui(pi)d*(p 2 ), (23) 

where i,j — 1,...,6 label up- and down-type squark mass eigenstates, h a b helicities, and p a ,b,i,2 four-momenta, 
proceeds from an equally charged quark-antiquark initial state through the tree-level Feynman diagrams shown in 
Fig. [U The corresponding cross section can be written in a compact way as 
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FIG. 1: Tree-level Feynman diagrams for the production of charged squark-antisquark pairs in quark-antiquark collisions. 
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which combine coupling constants and Dirac traces of the squared and interference diagrams. In cMFV, superpartners 
of heavy flavours can only be produced through the purely left-handed s-channel W-exchange, since the i-channel 
diagrams are supp ressed by the small bottom and negligible top quark densities in the proton, and one recovers the 
result in Ref. [20|. In NMFV, c-channel exchanges can, however, contribute to heavy-flavour final state production 
from light-flavour initial states and even become dominant, due to the strong gluino coupling. 

Neutral squark-antisquark pair production proceeds either from equally neutral quark-antiquark initial states 

q(ha,Pa)q{h b ,p b ) -> qi(pi)q*(p2), (26) 

through the five different gauge-boson/gaugino exchanges shown in Fig. [2] (top) or from gluon-gluon initial states 

g(h a ,p a )g(h b ,p b ) -> qi(pi)q*(p2) (27) 



through the purely strong couplings shown in Fig. [2] (bottom). The differential cross section for quark-antiquark 
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FIG. 2: Tree-level Feynman diagrams for the production of neutral squark-antisquark pairs in quark-antiquark (top) and 
gluon-gluon collisions (bottom). 
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FIG. 3: Tree-level Feynman diagrams for the production of one down-type squark (qi) and one up-type squark (q'j) in the 
collision of an up- type quark (q) and a down- type quark (g'). 
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since only very few interferences (those between strong and electroweak channels of the same propagator type) are 
eliminated due to colour conservation. On the other hand, the gluon-initiated cross section 
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involves only the strong coupling constant and is thus quite compact. In the case of cMFV, but diagonal or non- 
diagonal squark hclicity, our results agree with those in Ref. [26[ . Diagonal production of identical squark-antisquark 
mass eigenstates is, of course, dominated by the strong quark-antiquark and gluon-gluon channels. Their relative im- 
portance depends on the partonic luminosity and thus on the type and energy of the hadron collider under considera- 
tion. Non-diagonal production of squarks of different helicity or flavour involves only electroweak and gluino-mediated 
quark-antiquark scattering, and the relative importance of these processes depends largely on the gluino mass. 



C. Squark Pair Production 



While squark-antisquark pairs are readily produced in pp collisions, e.g. at the Tevatron, from valence quarks and 
antiquarks, pp colliders have a larger quark-quark luminosity and will thus more easily lead to squark pair production. 
The production of one down-type and one up-type squark 

q{h a ,p a )q'(h b ,p b ) -> di(pi) Uj{jp 2 ), (31) 

in the collision of an up-type quark q and a down-type quark q' proceeds through the i-channel chargino or it-channel 
neutralino and gluino exchanges shown in Fig. [3] The corresponding cross section 
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FIG. 4: Tree-level Feynman diagrams for the production of two up-type or down-type squarks. 
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(33) 



where the neutralino-gluino interference term is absent due to colour conservation. The cross section for the charge- 
conjugate production of antisquarks from antiquarks can be obtained from the equations above by replacing h a ^ b — > 
— h a b . Heavy-flavour final states are completely absent in cMFV due to the negligible top quark and small bottom 
quark densities in the proton and can thus only be obtained in NMFV. 

The Feynman diagrams for pair production of two up- or down-type squarks 



q(h a ,Pa)q'(h b ,p b ) -> q l (pi)q j {p2) 



(34) 



are shown in Fig. |U In NMFV, neutralino and gluino exchanges can lead to identical squark flavours for different 
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quark initial states, so that both t- and u-channels contribute and may interfere. The cross section 
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depends therefore on the form factors 



[AAT] 
[MU] 
[MTU] 

\GT\rr 

[GU]„ 

[GTU]rr 

WGA] k „ 
WGB] k 



kl 

mn 



kl 
mn 



kl 
mn 



it a 2 


r 


x 2 w (1 - x w 


2 s 2 


ira 2 


r 


x w (! _ x w 


2 S 2 U 


2ira 2 





n* nm* nn 



Zx 2 w (1 - x w ) 2 s 2 
2 ir a 2 

qjq'g <iiQ§ 



Re 



nm* nn* nn nm 
Qiixl Qjl'xl gil'x" Qjix" 



9s 2 

2wa 2 s 
9 s 2 

-in a 2 



nm nn 



Re 



27 s 2 
8 7raa 



9 s 2 xw (1 — xw) 



8 7r aa s 



9 s 2 xw (1 — xw) 



Re 



Re 



(ut - to? to|. ) (1 - <W) + "^^o "i^o « 5( 
to? m?. ) (1 — <W) + "^x!! m x? s ^ 
^ut - to? (1 - S rnn ) + to^o to x o s 5, 

(u t-m^.m-^ (1 - 5 mn ) + to? s <5 m 

f - TO? TO?J (1 - (5„ m ) + TO? S S m 

(ut- m| to?.) (1 - 5 mn ) + to? s 6- 

(lit- TO? TO? ) (1 - 5 mn ) + TO x o TOg S <5-, 

^ui - TO? TO? ) (1 - S mn ) + m^o m g S 5; 



n rn nn nm* nn* 
U qiqg ^qjq'g^qiq'g U q~jqg 



nn* nm* nm* nn* 

U Qji'x k u q,qx° k ^qiq'S^qjqS 



nn* nm* nn* nm* 

u q t q'xl u qjqx° k 1 <(;■/';/ 1 ■/ ■,"/ 



(36) 



Gluinos will dominate over neutralino exchanges due to their strong coupling, and the two will only interfere in the 
mixed t- and u-channels due to colour conservation. At the LHC, up-type squark pair production should dominate 
over mixed up-/down-type squark production and down- type squark pair production, since the proton contains two 
valence up-quarks and only one valence down-quark. As before, the charge-conjugate production of antisquark pairs 
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FIG. 5: Tree-level Feynman diagrams for the associated production of squarks and gauginos. 



is obtained by making the replacement h a ,b —* —h a ,b- If we neglect electroweak contributions as well as squark flavour 
and helicity mixing and sum over left- and right-handed squark states, our results agree with those of Ref. [27| . 



D. Associated Production of Squarks and Gauginos 



The associated production of squarks and neutralinos or charginos 

q(ha,Pa)g(hb,Pb) -> Xj(pi)Qi(P2) 



(37) 



is a semi- weak process that originates from quark-gluon initial states and has both an s-channel quark and a i-channel 
squark contribution. They involve both a quark-squark-gaugino vertex that can in general be flavour violating. The 
corresponding Feynman diagrams can be seen in Fig. [5) The squark-gaugino cross section 
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(38) 



where 



6xw(l — %w) for neutralinos and — 12xw for charginos, is sufficiently compact to be written without 



the definition of form factors. Note that the i-channel diagram involves the coupling of the gluon to scalars and 
does thus not depend on its helicity h^. The cross section of the charge-conjugate process can be obtained by taking 
h a — > — h a . Third-generation squarks can only be produced in NMFV, preferably through a light (valence) quark in 
the s-channel. For non- mixing squarks and gauginos, we agree again with the results of Ref. [271 ]. 



E. Gaugino Pair Production 

Finally, we consider the purely electroweak production of gaugino pairs 

q(h a ,Pa)q'(h b ,p b ) -> Xi(Pl)Xj(P2) (39) 

from quark-antiquark initial states, where flavour violation can occur via the quark-squark-gaugino vertices in the t- 
and u-channels (see Fig. However, if it were not for different parton density weights, summation over complete 
squark multiplet exchanges would make these channels insensitive to the exchanged squark flavour. Furthermore there 
are no final state squarks that could be experimentally tagged. The cross section can be expressed generically as 
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FIG. 6: Tree-level Feynman diagrams for the production of gaugino pairs. 
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(40) 



i.e. in terms of generalized charges. For \i -production, these charges are given by 
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(41) 



Note that there is no interference between t- and it-channel diagrams due to (electromagnetic) charge conservation. 
The cross section for chargino-pair production in e + e~-collisions can be deduced by setting e q — > e; = — 1, L qq 'z ~ * 
^eeZ — (2 T/ 3 — 2 e; xy) and R qq 'z ~ * ReeZ = —2 e/ xvk- Neglecting all Yukawa couplings, we can then reproduce the 
calculations of Ref . [28| ■ 
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The charges of the chargino-neutralino associated production are given by 
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(42) 



The charge-conjugate process is again obtained by making the replacement h a ,b — > — h a ^ in Eq. (|40| . In the case of 
non- mixing squarks with neglected Yukawa couplings, we agree with the results of Ref. 22], provided we correct a 
sign in their Eq. (2) as described in Ref. [29| . 

Finally, the charges for the neutralino pair production are given by 
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(43) 



which agrees with the results of Ref. [30| in the case of non-mixing squarks. 
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FIG. 7: Tree-level Feynman diagrams for squark decays into gauginos and quarks (top) and into electroweak gauge bosons and 
lighter squarks (bottom). 



F. Squark Decays 



We turn now from SUSY particle production to decay processes and show in Fig. the possible decays of squarks 
into gauginos and quarks (top) as well as into electroweak gauge bosons and lighter squarks (bottom). Both processes 
can in general induce flavour violation. The decay widths of the former are given by 
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In cMFV, our results agree with those of Ref. [3l[ . 
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G. Gluino Decays 



Heavy gluinos can decay strongly into squarks and quarks as shown in Fig. [5] The corresponding decay width 
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FIG. 8: Tree-level Feynman diagram for gluino decays into squarks and quarks 






FIG. 9: Tree-level Feynman diagrams for gaugino decays into squarks and quarks (left) and into lighter gauginos and electroweak 
gauge bosons (centre and right). 
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H. Gaugino Decays 

Heavier gauginos can decay into squarks and quarks as shown in Fig. [9] (left) or into lighter gauginos and electroweak 
gauge bosons (Fig. [9] centre and right). The analytical decay widths are 
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(56) 



for neutralinos, respectively. Chargino decays into a slepton and a neutrino (lepton and sneutrino) can be deduced 
from the previous equations by taking the proper limits, i.e. by removing colour factors and up-type masses in the 
coupling definitions. Our results agree then with those of Ref. [32| in the limit of non- mixing sneutrinos. Note that 
the same simplifications also permit a verification of our results for squark decays into a gaugino and a quark in Eqs. 
(|44|) and (|45|) when compared to their leptonic counterparts in Ref. [32| ■ 



IV. EXPERIMENTAL CONSTRAINTS, SCANS AND BENCHMARKS IN NMFV SUSY 



In the absence of experimental evidence for supersymmetry, a large variety of data can be used to constrain the 
MSSM parameter space. For example, sparticle mass limits can be obtained from searches of charginos (m~± > 85 

Xi 

GeV for heavier sneutrinos at LEP2), neutralinos (m^o > 59 GeV in minimal supergravity (mSUGRA) from the 
combination of LEP2 results), gluinos (rrig > 195 GeV from CDF), stops (rnj > 95... 96 GeV for neutral- or 
charged-current decays from the combination of LEP2 results), and other squarks (rriq > 300 GeV for gluinos of equal 
mass from CDF) at colliders [33j. Note that all of these limits have been obtained assuming minimal flavour violation. 
For non-minimal flavour violation, rather strong constraints can be obtained from low-energy, electroweak precision, 
and cosmological observables. These are discussed in the next subsection, followed by several scans for experimentally 
allowed/favoured regions of the constrained MSSM parameter space and the definition of four NMFV benchmark 
points/slopes. Finally, we exhibit the corresponding chirality and flavour decomposition of the various squark mass 
eigenstates. 



A. Low-Energy, Electroweak Precision, and Cosmological Constraints 

In a rather complete analysis of FCNC constraints more than ten years ago upper limits from the neutral 
kaon sector (on ArriK, £, e'/e), on B- (Arris) and D-meson oscillations (Amu), various rare decays (BR(& — > sy), 
BR(/i — > 67), BR(r ~> e*y), and BR(r — > fJ/y)), and electric dipole moments (d n and d e ) were used to impose 
constraints on non-minimal flavour mixing in the squark and slepton sectors. The limit obtained for the absolute 
value in the left-handed, down-type squark sector was rather weak (\\f L \ < 4.4... 26 for varying gluino-to-squark 
mass ratio), while the limits for the mixed/right-handed, imaginary or sleptonic parts were already several orders of 
magnitude smaller. In the meantime, many of the experimental bounds have been improved or absolute values for the 
observables have been determined, so that an updated analysis could be performed [34| . The results for the down-type 
squark sector are cited in Tab. |TJ As can be seen and as has already been hinted at in the introduction, only mixing 
between second- and third-generation squarks can be substantial, and this only in the left-left or right-right chiral 
sectors, the latter being disfavoured by its scaling with the soft SUSY-breaking mass. Independent analyses focusing 
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TABLE I: The 95% probability bounds on |Afj= d, | obtained in Ref. JH- 
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on this particular sector, i.e. on BR(6 — > sj), BR(6 — > sfifj,), and AmB a , have been performed recently by two other 
groups ,35; .'!() with very similar results. 

In our own analysis, we take implicitly into account all of the previously mentioned constraints by restricting 
ourselves to the case of only one real NMFV parameter, A = \ S £ L — A|f L . Allowed regions for this parameter are then 
obtained by imposing explicitly a number of low-energy, electroweak precision, and cosmological constraints. We start 
by imposing the theoretically robust inclusive branching ratio 

BR(6^ si) = (3.55 ±0.26) x 1CT 4 , (57) 

obtained from the combined measurements of BaBar, Belle, and CLEO jsjj, at the 2cr-level on the two-loop QCD/one- 
loop SUSY calculation [36|, [38|, which affects directly the allowed squark mixing between the second and third 
generation. 

A second important consequence of NMFV in the MSSM is the generation of large splittings between squark- mass 
eigenvalues. The splitting within isospin doublets influences the Z- and lU-boson self-energies at zero-momentum 
£ z,w(0) in the electroweak p-parameter 

Ap = S z (0)/M|-S w (0)/M^ (58) 

and consequently the lU-boson mass My/ and the squared sine of the weak mixing angle sin 2 8\y- The latest combined 
fits of the Z-boson mass, width, pole asymmetry, IF-boson and top-quark mass constrain new physics contributions 
to T = -0.13 ± 0.11 H| or 

Ap = -oT = 0.00102 ± 0.00086, (59) 

where we have used a(Mz) = 1/127.918. This value is then imposed at the 2cr-level on the one-loop NMFV and 
two-loop cMFV SUSY calculation [HI]. 

A third observable sensitive to SUSY loop-contributions is the anomalous magnetic moment a M = (g^ — 2)/2 of the 
muon, for which recent BNL data and the SM prediction disagree by [H| 

Aet M = (22 ± 10) x 10 -10 . (60) 

In our calculation, we take into account the SM and MSSM contributions up to two loops [I(J EH and require them 
to agree with the region above within two standard deviations. 

For cosmological reasons, i.e. in order to have a suitable candidate for non-baryonic cold dark matter [42j], we require 
the lightest SUSY particle (LSP) to be stable, electrically neutral, and a colour singlet. The dark matter relic density 
is then calculated using a modified version of DarkSUSY 4.1 [43[, that takes into account the six-dimensional squark 
helicity and flavour mixing, and constrained to the region 

0.094 < n CDM h 2 < 0.136 (61) 

at 95% (2a) confidence level. This limit has recently been obtained from the three-year data of the WMAP satellite, 
combined with the SDSS and SNLS survey and Baryon Acoustic Oscillation data and interpreted within an eleven- 
parameter inflationary model [4^ |. which is more general than the usual six-parameter "vanilla" concordance model 
of cosmology. Note that this range is well compatible with the older, independently obtained range of 0.094 < 
ttcDMh 2 < 0.129 [H. 

B. Scans of the Constrained NMFV MSSM Parameter Space 

The above experimental limits are now imposed on the constrained MSSM (cMSSM), or minimal supergravity 
(mSUGRA), model with five free parameters mo, m\j2, tan/3, Aq, and sgn(/x) at the grand unification scale. Since 
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FIG. 10: The (m , m 1/2 )-planes for tan/3 = 10, A = GeV, fi < 0, and A = 0, 0.03, 0.05 and 0.1. We show WMAP (black) 
favoured as well as b — > sy (blue) and charged LSP (beige) excluded regions of mSUGRA parameter space in minimal (A = 0) 
and non-minimal (A > 0) flavour violation. 

our scans of the cMSSM parameter space in the common scalar mass mo and the common fermion mass m.1/2 depend 
very little on the trilinear coupling Aq, we set it to zero in the following. Furthermore, we fix a small (10), intermediate 
(30), and large (50) value for the ratio of the Higgs vacuum expectation values tan/3. The impact of the sign of the 
off-diagonal Higgs mass parameter fi is investigated for tan/3 = 10 only, before we set it to fi > for tan/3 = 30 and 
50 (see below). 

With these boundary conditions at the grand unification scale, we solve the renormalization group equations numer- 
ically to two-loop order using the computer program SPheno 2.2.3 [!(| and compute the soft SUSY-breaking masses 
at the electroweak scale with the complete one-loop formulas, supplemented by two-loop contributions in the case of 
the neutral Higgs bosons and the /i-parameter. At this point we generalize the squark mass matrices as described 
in Sec. |TT]in order to account for flavour mixing in the left-chiral sector of the second- and third-generation squarks, 
diagonalize these mass matrices, and compute the low-energy, electroweak precision, and cosmological observables 
with the computer programs FeynHiggs 2.5.1 [47| an d DarkSUSY 4.1 [43l |. 

For the masses and widths of the electroweak gauge bosons and the mass of the top quark, we use the current 
values of m z = 91.1876 GeV, m w = 80.403 GeV, m t = 174.2 GeV, T z = 2.4952 GeV, and T w = 2.141 GeV. The 
CKM-matrix elements are computed using the parameterization 

(C12C13 S12C13 si S e- %s \ 

-S12C23 - c 12 s 2 3Si3e tS C12C23 - s 12 s 2 3Si3e lS s 23 ci3 , (62) 
S12S23 - Ci 2 C2 3 Sl3e J ' 5 -C12S23 - Si 2 C23Sl 3 e 1 ' 5 C23C13 / 

where = sin Oij and cy = cos Oij relate to the mixing of two specific generations i and j and 8 is the SMC P- violating 
complex phase. The numerical values are given by 

S12 = 0.2243, s 23 = 0.04 1 3, s 13 = 0.0037, and 8 = 1.05. (63) 

The squared sine of the electroweak mixing angle sin 2 0w = 1 — rr ^w I m % anc ^ ^ ne electromagnetic fine structure 
constant a — \[2G fin^ sin 2 9w/^ are calculated in the improved Born approximation using the world average value 
of Gf = 1.16637 ■ 10~ 5 GeV~ 2 for Fermi's coupling constant (33[. 

Typical scans of the cMSSM parameter space in mo and m^ 2 with a relatively small value of tan (3 = 10 and Aq = 
are shown in Figs. [T0l and [TT1 for /.i < and /1 > 0, respectively. All experimental limits described in Sec. II V Al arc 
imposed at the 2<r-level. The b — > sj excluded region depends strongly on flavour mixing, while the regions favoured 
by 3/x — 2 and the dark matter relic density are quite insensitive to variations of the A-parameter. Ap constrains the 
parameter space only for heavy scalar masses mo > 2000 GeV and heavy gaugino masses mi/2 > 1500 GeV, so that 
the corresponding excluded regions are not shown here. 

The dominant SUSY effects in the calculation of the anomalous magnetic moment of the muon come from induced 
quantum loops of a gaugino and a slepton. Squarks contribute only at the two-loop level. This reduces the dependence 
on flavour violation in the squark sector considerably. Furthermore, the region /x < is disfavoured in all SUSY models, 
since the one-loop SUSY contributions are approximatively given by [48[ 

a susY, i-!oo P ^ 13 x 10 -io tan/3 sgn( M ), (64) 

if all SUSY particles (the relevant ones are the smuon, sneutralino, chargino, and neutralino) have a common mass 
-^susy- Negative values of /i would then increase, not decrease, the disagreement between the experimental measure- 
ments and the theoretical SM value of a^. Furthermore, the measured b — > sj branching ratio excludes virtually all 
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FIG. 11: The (m , m 1/2 )-planes for tan/3 = 10, A = GeV, /x > 0, and A = 0, 0.03, 0.05 and 0.1. We show a M (grey) and 
WMAP (black) favoured as well as b — > s'y (blue) and charged LSP (beige) excluded regions of mSUGRA parameter space in 
minimal (A = 0) and non-minimal (A > 0) flavour violation. 
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FIG. 12: The (m ,m 1/2 ) planes for tan/3 = 30, A = GeV, /Lt > 0, and A = 0, 0.03, 0.05 and 0.1. We show a M (grey) and 
WMAP (black) favoured as well as 6 — » S7 (blue) and charged LSP (beige) excluded regions of mSUGRA parameter space in 
minimal (A = 0) and non-minimal (A > 0) flavour violation. 



of the region favoured by the dark matter relic density, except for very high scalar SUSY masses. We therefore do 
not consider negative values of [i in the rest of this work. As stated above, we have also checked that the shape of the 
different regions depends extremely weakly on the trilinear coupling Aq. 

In Figs. [T!?l and [TUI we show the (mo, m 1 / 2 )-planes for larger tan/3, namely tan/3 = 30 and tan/3 = 50, and for 
/i > 0. The regions which are favoured both by the anomalous magnetic moment of the muon and by the cold dark 
matter relic density, and which are not excluded by the b — > S7 measurements, are stringently constrained and do not 
allow for large flavour violation. 
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FIG. 13: The (m ,m 1/2 ) planes for tan/3 = 50, A = GeV, /i > 0, and A = 0, 0.03, 0.05 and 0.1. We show a M (grey) and 
WMAP (black) favoured as well as 6 — » s'y (blue) and charged LSP (beige) excluded regions of mSUGRA parameter space in 
minimal (A = 0) and non-minimal (A > 0) flavour violation. 
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TABLE II: Benchmark points allowing for flavour violation among the second and third generations for Ao = 0, (X > 0, and 



three different values of tan/3. For comparison we also show the nearest pre-WMAP SPS 
[5lj benchmark points and indicate the relevant cosmological regions. 



50] and post-WMAP BDEGOP 



m [GeV] m 1/2 [GeV] A [GeV] tan/j sgn(/i) 
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C. (c)MFV and NMFV Benchmark Points and Slopes 

Restricting ourselves to non-negative values of /i, we now inspect the (too, TO!/ 2 )-planes in Figs. fTTTll3l for cMSSM 
scenarios that 

• are allowed/favoured by low-energy, electroweak precision, and cosmological constraints, 

• permit non-minimal flavour violation among lcft-chiral squarks of the second and third generation up to A < 0.1, 

• and are at the same time collider- friendly, i.e. have relatively low values of mo and m 1 / 2 . 

Our choices are presented in Tab.Hl] together with the nearest pre-WMAP Snowmass Points (and Slopes, SPS) (49l.l50l] 
and the nearest post-WMAP scenarios proposed in Ref. [5l|. We also indicate the relevant cosmological region for 
each point and attach a model line (slope) to it, given by 



A : 


180 GeV 


< 




< 


250 GeV , 


to = 


- 1936 GeV + 


12.9to 1/2 


B : 


400 GeV 


< 


mi/2 


< 


900 GeV , 


to = 


4.93 GeV + 


0.229 toi/2 


C : 


500 GeV 


< 


m 1 / 2 


< 


700 GeV , 


to = 


54 GeV + 


0.297 mi/2 


D : 


575 GeV 


< 


mi/2 


< 


725 GeV , 


to = 


600 GeV. 





These slopes trace the allowed/favoured regions from lower to higher masses and can, of course, also be used in cMFV 

scenarios with A = 0. We have verified that in the case of MFV Q the hierarchy A'l ^> A'^ RL A RR and the 

equality of A£ £ = A^ L are still reasonably well fulfilled numerically with the values of \ S £ L ~ X^ L ranging from zero 
to5xl0~ 3 ...lxl0~ 2 for our four typical benchmark points. 

Starting with Fig. [TT1 and tan/? = 10, the bulk region of equally low scalar and fermion masses is all but excluded 
by the b — > sj branching ratio. This leaves as a favoured region first the so-called focus point region of low fermion 
masses m 1 / 2 , where the lightest neutralinos are relatively heavy, have a significant Higgsino component, and annihilate 
dominantly into pairs of electroweak gauge bosons. Our benchmark point A lies in this region, albeit at smaller masses 
than SPS 2 (to = 1450 GeV, to 1/2 = 300 GeV) and BDEGOP E' (to = 1530 GeV, m 1/2 = 300 GeV), which lie 
outside the region favoured by (grey-shaded) and lead to collider-unfriendly squark and gaugino masses. 

The second favoured region for small tan f3 is the co-annihilation branch of low scalar masses too, where the lighter 
tau-slepton mass eigenstate is not much heavier than the lightest neutralino and the two have a considerable co- 
annihilation cross section. This is where we have chosen our benchmark point B, which differs from the points SPS 
3 (to = 90 GeV, to 1/2 = 400 GeV) and BDEGOP C' (m = 85 GeV, to 1/2 = 400 GeV) only very little in the scalar 
mass. This minor difference may be traced to the fact that we use DarkSUSY 4.1 [43| instead of the private dark 
matter program SSARD of Ref. [5l| . 

At the larger value of tan [3 = 30 in Fig. 1121 only the co-annihilation region survives the constraints coming from 
b — > S7 decays. Here we choose our point C, which has slightly higher masses than both SPS lb (too = 200 GeV, 
m i/2 — 400 GeV) and BDEGOP F (too = 175 GeV, mi/2 = 350 GeV), due to the ever more stringent constraints 
from the above-mentioned rare -B-decay. 

For the very large value of tan/3 = 50 in Fig. 1131 the bulk region reappears at relatively heavy scalar and fermion 
masses. Here, the couplings of the heavier scalar and pseudo-scalar Higgses H° and ^4° to bottom quarks and tau- 
lcptons and the charged-Higgs coupling to top-bottom pairs are significantly enhanced, resulting e.g. in increased 
dark matter annihilation cross sections through s-channel Higgs-exchange into bottom-quark final states. So as tan (3 
increases further, the so-called Higgs- funnel region eventually makes its appearance on the diagonal of large scalar and 
fermion masses. We choose our point D in the concentrated (bulky) region favoured by cosmology and a M at masses, 
that are slightly higher than those of SPS 4 (to = 400 GeV, to 1/2 = 300 GeV) and BDEGOP L' (to = 300 GeV, 
m i/2 = 450 GeV). We do so in order to escape again from the constraints of the b — > sj decay, which are stronger 
today than they were a few years ago. In this scenario, squarks and gluinos are very heavy with masses above 1 TeV. 



21 





0.2 0.4 0.6 O.S 1 0.2 0.4 0.6 OS 1 0.2 0.4 0.6 OS 



FIG. 14: Dependence of the precision variables BR(fe — > sj), Ap, and the cold dark matter relic density QcDnih 2 (top) as well 
as of the lightest SUSY particle, up- and down-type squark masses (bottom) on the NMFV parameter A in our benchmark 
scenario A. The experimentally allowed ranges (within 2a) are indicated by horizontal dashed lines. 



D. Dependence of Precision Observables and Squark-Mass Eigenvalues on Flavour Violation 

Let us now turn to the dependence of the precision variables discussed in Sec. lIV Al on the flavour violating parameter 
A in the four benchmark scenarios defined in Sec. IIV CI As already mentioned, we expect the leptonic observable 
to depend weakly (at two loops only) on the squark sector, and this is confirmed by our numerical analysis. We find 
constant values of 6, 14, 16, and 13xl0 -10 for the benchmarks A, B, C, and D, all of which lie well within 2a (the 
latter three even within la) of the experimentally favoured range (22 ± 10) x 1CP 10 . 

The electroweak precision observable Ap is shown first in Figs. fT4lfT71 for the four benchmark scenarios A, B, C, and 
D. On our logarithmic scale, only the experimental upper bound of the 2cr-range is visible as a dashed line. While the 
self-energy diagrams of the electroweak gauge bosons depend obviously strongly on the helicities, flavours, and mass 
eigenvalues of the squarks in the loop, the SUSY masses in our scenarios are sufficiently small and the experimental 
error is still sufficiently large to allow for relatively large values of A < 0.57, 0.52, 0.38, and 0.32 for the benchmark 
points A, B, C, and D, respectively. As mentioned above, Ap conversely constrains SUSY models in cMFV only for 
masses above 2000 GeV for too an d 1500 GeV for m 1 / 2 - 

The next diagram in Figs. 11411 1 71 shows the dependence of the most stringent low-energy constraint, coming from 
the good agreement between the measured b — > S7 branching ratio and the two-loop SM prediction, on the NMFV 
parameter A. The dashed lines of the 2<r-bands exhibit two allowed regions, one close to A = (vertical green line) 
and a second one around A ~ 0.57, 0.75, 0.62, and 0.57, respectively. As is well-known, the latter are, however, 
disfavoured by b — > sp + p~ data constraining the sign of the b — > amplitude to be the same as in the SM [H^]. We 
will therefore limit ourselves later to the regions A < 0.05 (points A, C, and D) and A < 0.1 (point B) in the vicinity 
of (c)MFV (see also Tab.H|. 

The 95% confidence-level (or 2a) region for the cold dark matter density was given in absolute values in Ref. [3| 
and is shown as a dashed band in the upper right part of Figs. [MlflTl However, only the lower bound (0.094) is of 
relevance, as the relic density falls with increasing A. This is not so pronounced in our model B as in our model A, 
where squark masses are light and the lightest ncutralino has a sizable Higgsino-component, so that squark exchanges 
contribute significantly to the annihilation cross sections. For models C and D there is little sensitivity of SlcDMh 2 
(except at very large A < 1), as the squark masses are generally larger. 




FIG. 16: Same as Fig. 1141 for our benchmark scenario C. 
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FIG. 17: Same as Fig. [T4]for our benchmark scenario D. 



The rapid fall-off of the relic density for very large A < 1 can be understood by looking at the resulting lightest up- 
and down-type squark mass eigenvalues in the lower left part of Figs. [MlfTTl For maximal flavour violation, the off- 
diagonal squark mass matrix elements are of similar size as the diagonal ones, leading to one squark mass eigenvalue 
that approaches and finally falls below the lightest neutralino (dark matter) mass. Light squark propagators and 
co-annihilation processes thus lead to a rapidly falling dark matter relic density and finally to cosmologically excluded 
NMFV SUSY models, since the LSP must be electrically neutral and a colour singlet. 

An interesting phenomenon of level reordering between neighbouring states can be observed in the lower central 
diagrams of Figs. [FjlfTTl for the two lowest mass eigenvalues of up-type squarks. The squark mass eigenstates are, by 
definition, labeled in ascending order with the mass eigenvalues, so that u\ represents the lightest, «2 the second- 
lightest, and uq the heaviest up-type squark. As A and the off-diagonal entries in the mass matrix increase, the splitting 
between the lightest and highest mass eigenvalues naturally increases, whereas the intermediate squark masses (of 
^3,4,5) are practically degenerate and insensitive to A. These remarks also hold for the down-type squark masses 
shown in the lower right diagrams of Figs. I14II17I However, for up-type squarks it is first the second-lowest mass that 
decreases up to intermediate values of A = 0.2. ..0.5, whereas the lowest mass is constant, and only at this point the 
second-lowest mass becomes constant and takes approximately the value of the until here lowest squark mass, whereas 
the lowest squark mass starts to decrease further with A. These "avoided crossings" are a common phenomenon for 
Hermitian matrices and reminiscent of meta-stable systems in quantum mechanics. At the point where the two 
levels should cross, the corresponding squark eigenstates mix and change character, as will be explained in the next 
subsection. For scenario C (Fig. I16|) , the phenomenon occurs even a second time with an additional avoided crossing 
between the states U2 and U3 at A ~ 0.05. For scenario B (Fig. [T5| . this takes place at A ~ 0.1, and there is even 
another crossing at A ~ 0.02. For down-type squarks, the level-reordering phenomenon is not so pronounced. 



E. Chirality and Flavour Decomposition of Squark Mass Eigenstates 

In NMFV, squarks will not only exhibit the traditional mixing of left- and right-handed helicities of third-generation 
flavour eigenstates, but will in addition exhibit generational mixing. As discussed before, we restrict ourselves here to 
the simultaneous mixing of left-handed second- and third-generation up- and down-type squarks. For our benchmark 
scenario A, the helicity and flavour decomposition of the six up-type (left) and down-type (right) squark mass eigen- 




FIG. 19: Same as Fig.lTHlfor A e [0;0.1]. 




FIG. 21: Same as Fig.l2Qlfor A G [0;0.1]. 




FIG. 23: Same as Fig.l22lfor A G [0;0.1]. 




FIG. 25: Same as Fig. HU for A G [0;0.1]. 
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states is shown in Fig. [18] for the full range of the parameter A € [0; 1] and in Fig. [19] for the experimentally favoured 
range in the vicinity of (c)MFV, A G [0; 0.1]. First-generation and right-handed second-generation squarks remain, of 
course, helicity- and flavour-diagonal, 

"3 = c R , d 2 = s R , 
"4 = u R , d 3 = d R , 

u 5 = Ul , cfe = d~L, (66) 

with the left-handed and first-generation squarks being slightly heavier due their weak isospin coupling (see Eqs. ([5])- 
([?])) and different renormalization-group running effects. Production of these states will benefit from t- and it-channel 
contributions of first- and second-generation quarks with enhanced parton densities in the external hadrons, but they 
will not be identified easily with heavy-flavour tagging and are of little interest for our study of flavour violation. 
The lightest up-type squark u\ remains the traditional mixture of left- and right-handed stops over a large region of 
A < 0.4, but it shows at this point an interesting flavour transition, which is in fact expected from the level reordering 
phenomenon discussed in the lower central plot of Fig. 1141 The transition happens, however, above the experimental 
limit of A < 0.1. Below this limit, it is the states U2, uq, d±, and in particular gLi and de that show, in addition to 
helicity mixing, the most interesting and smooth variation of second- and third-generation flavour content (see Fig. 
[T9]) . Note that at very low A ~ 0.002 the states d^ and sl rapidly switch levels. This numerically small change was 
not visible on the linear scale in Fig. [TJ] 

For the benchmark point B, whose helicity and flavour decomposition is shown in Fig. I20| level reordering occurs 
at A ~ 0.4 for the intermediate-mass up-type squarks, 



"3,4 


= c R 


, d~2 


= Sr. 


"4,3 


= u R 


, d~3 


= d R 


Us 


= u L 


, d~5 


= d L 



(67) 

whereas the ordering of down- type squarks is very similar to scenario A. Close inspection of Fig. [5l] shows, however, 
that also d R and s R switch levels at low values of A ~ 0.02. At A ~ 0.01, in addition s R and switch levels, and at 
A ~ 0.002 it is the states ul and cl- The lightest up-type squark is again nothing but a mix of left- and right-handed 
stops up to A < 0.4. Phenomenologically smooth transitions below A < 0.1 involving taggable third-generation squarks 
are observed for -S4, uq, d%, and d$. 

The helicity and flavour decomposition for our scenario D, shown in Fig. 1241 is rather similar to the one in scenario 
A, i.e. 

"3 = Cr , d 3 = s R , 
u A = u R , J 4 = d R , 

u 5 = u L , d 5 — d L (68) 

are exactly the same in the up-squark sector, and only the mixed down- type state d^ is now lighter and becomes d 2 . 
The lightest up-type squark, Gi, is again mostly a mix of left- and right-handed top squarks up to A ~ 0.4, where 
the level reordering and generation mixing occurs (see lower central part of Fig. 1 1 7[) . At the experimentally favoured 
lower values of A < 0.1, the states U2, tie, di, d 2 , and d§ exhibit some smooth variations, shown in detail in Fig. [25] 
albeit to a lesser extent than in scenario A. At very low A ~ 0.004, it is now the up-type squarks ul and cl that 
rapidly switch levels. This numerically small change was again not visible on a linear scale (see Fig. H7|) . 
For our scenario C, shown in Fig. [5H the assignment of the intermediate states 

"3 = c R , d 3 = SR, 
"4 = ur , d A — d R , 

u 5 = u L , d 5 = d L (69) 

is the same as for scenario D above A > 0.1. Just below, u R and c R as well as d R and s R rapidly switch levels, whereas 
ul and cl switch levels at very low A ~ 0.002. These changes were already visible upon close inspection of the lower 
central and right plots in Fig. 1161 On the other hand, the lightest squarks u\ and d\ only acquire significant flavour 
admixtures at relatively large A ~ 0.2. ..0.4, whereas they are mostly superpositions of left- and right-handed stops 
and sbottoms in the experimentally favourable range of A < 0.1 shown in Fig. [23] Here, the heaviest states uq and d§ 
show already smooth admixtures of third-generation squarks as it was the case for the scenarios A and D discussed 
above. The most interesting states are, however, u 2 , U4, d 2 , and d±, respectively, since they represent combinations 
of up to four different helicity and flavour states and have a significant, taggable third-generation flavour content. 
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V. NUMERICAL PREDICTIONS FOR NMFV SUSY PARTICLE PRODUCTION AT THE LHC 

In this section, we present numerical predictions for the production cross sections of squark-antisquark pairs, squark 
pairs, the associated production of squarks and gauginos, and gaugino pairs in NMFV SUSY at the CERN LHC, i.e. 
for pp-collisions at y/S = 14 TeV centre-of-mass energy. Thanks to the QCD factorization theorem, total unpolarized 
hadronic cross sections 



can be calculated by convolving the relevant partonic cross sections da/dt, computed in Sec. lIIIi with universal parton 
densities f a /A an d fb/B of partons a,b in the hadrons A,B, which depend on the longitudinal momentum fractions 
of the two partons x a ^ = y/re ±y and on the unphysical factorization scales M a fi- For consistency with our leading 
order (LO) QCD calculation in the collinear approximation, where all squared quark masses (except for the top-quark 
mass) to 2 <§C s, we employ the LO set of the latest CTEQ6 global parton density fit [Hj], which includes nj — 5 "light" 
(including the bottom) quark flavours and the gluon, but no top-quark density. Whenever it occurs, i.e. for gluon 
initial states and gluon or gluino exchanges, the strong coupling constant a s ((iii) is calculated with the corresponding 
LO value of A^q - =165 MeV. We identify the renormalization scale /jr with the factorization scales M a = M& and 
set the scales to the average mass of the final state SUSY particles i and j, m = (m, + mj)/2. 

The numerical cross sections for charged squark-antisquark and squark-squark production, neutral up- and down- 
type squark-antisquark and squark-squark pair production, associated production of squarks with charginos and 
neutralinos, and gaugino pair production are shown in Fig. [2)] for our benchmark scenario A, in Fig. [23 for scenario 
B, in Fig. [221 for scenario C, and in Fig. [52] for scenario D. The magnitudes of the cross sections vary from the barely 
visible level of 10 -2 fb for weak production of heavy final states over the semi-strong production of average squarks 
and gauginos and quark-gluon initial states to large cross sections of 10 2 to 10 3 fb for the strong production of diagonal 
squark-(anti)squark pairs or weak production of very light gaugino pairs. Unfortunately, these processes, whose cross 
sections are largest (top right, center left, and lower right parts of Figs. |2"6"H2"9"]) . are practically insensitive to the flavour 
violation parameter A, as the strong gauge interaction is insensitive to quark flavours and gaugino pair production 
cross sections are summed over exchanged squark flavours. 

Some of the subleading, non-diagonal cross sections show, however, sharp transitions, in particular down-type 
squark-antisquark production at the benchmark point B (centre-left part of Fig. [27]), but also the other squark- 
antisquark and squark-squark production processes. At A = 0.02, the cross sections for d\d% and d^d% switch places. 
Since the concerned sstrange and sbottom mass differences are rather small, this is mainly due to the different strange 
and bottom quark densities in the proton. The cross section is mainly due to the exchange of strongly coupled gluinos 
despite their larger mass. At A = 0.035 the cross sections for d^d^ and did% increase sharply, since ds = dn can then 
be produced from down-type valence quarks. The cross section of the latter process increases with the strange squark 
content of d\. 

At the benchmark point C (Fig. [2"5]), sharp transitions occur between the U2/U4 and dijd^ states, which are pure 
charm/strange squarks below/above A = 0.035, for all types of charged and neutral squark-antisquark and squark- 
squark production and also squark-gaugino associated production. As a side-remark we note that an interesting 
perspective might be the exploitation of these i-channel contributions to second- and third-generation squark pro- 
duction for the determination of heavy-quark densities in the proton. This requires, of course, efficient experimental 
techniques for heavy-flavour tagging. 

Smooth transitions and semi-strong cross sections of about 1 fb are observed for the associated production of third- 
generation squarks with charginos (lower left diagrams) and neutralinos (lower centre diagrams) and in particular 
for the scenarios A and B. For benchmark point A (Fig. [26]) . the cross section for d^ production decreases with its 
strange squark content, while the bottom squark content increases at the same time. For benchmark point B (Fig. 
[27)) . the same (opposite) happens for de (di), while the cross sections for u§ increase/decrease with its charm/top 
squark content. Even in minimal flavour violation, the associated production of stops and charginos is a particularly 
useful channel for SUSY particle spectroscopy, as can be seen from the fact that cross sections vary over several orders 
of magnitude among our four benchmark points (see also Ref. (H)]). 

An illustrative summary of flavour violating hadroproduction cross section contributions for third-generation squarks 
and/or gauginos is presented in Tab. Hill together with the competing flavour-diagonal contributions. 
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FIG. 26: Cross sections for charged squark-antisquark (top left) and squark-squark (top centre) production, neutral up-type (top 
right) and down-type (centre left) squark-antisquark and squark-squark pair (centre and centre right) production, associated 
production of squarks with charginos (bottom left) and neutralinos (bottom centre), and gaugino pair production (bottom 
right) at the LHC in our benchmark scenario A. 



VI. CONCLUSIONS 



In conclusion, we have performed an extensive analysis of squark and gaugino hadroproduction and decays in non- 
minimal flavour violating supersymmetry. Within the super-CKM basis, we have taken into account the possible 
misalignment of quark and squark rotations and computed all squared helicity amplitudes for the production and the 
decay widths of squarks and gauginos in compact analytic form, verifying that our results agree with the literature in 
the case of non-mixing squarks whenever possible. Flavour violating effects have also been included in our analysis 
of dark matter (co-) annihilation processes. We have then analyzed the NMFV SUSY parameter space for regions 
allowed by low-energy, electroweak precision, and cosmological data and defined four new post-WMAP benchmark 
points and slopes equally valid in minimal and non-minimal flavour violating SUSY. We found that left-chiral mixing 
of second- and third-generation squarks is slightly stronger constrained than previously believed, mostly due to smaller 
experimental errors on the b — ► sj branching ratio and the cold dark matter relic density. For our four benchmark 
points, we have presented the dependence of squark mass eigenvalues and the flavour and helicity decomposition of 
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TABLE III: Dominant s-, t-, and it-channel contributions to the flavour violating hadroproduction of third-generation squarks 
and/or gauginos and the competing dominant flavour-diagonal contributions. 
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FIG. 28: Same as Fig. 1261 for our benchmark scenario C. 



the squark mass eigenstates on the flavour violating parameter A. We have computed numerically all production 
cross sections for the LHC and discussed in detail their dependence on flavour violation. A full experimental study 
including heavy-flavour tagging efficiencies, detector resolutions, and background processes would, of course, be very 
interesting in order to establish the experimental significance of NMFV. While the implementation of our analytical 
results in a general-purpose Monte Carlo generator should now be straight-forward, such a detailed experimental 
study represents a research project of its own [551 ] and is beyond the scope of the work presented here. 
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FIG. 29: Same as Fig. [26]for our benchmark scenario D. 



APPENDIX A: GAUGINO AND HIGGSINO MIXING 

The soft SUSY-brcaking terms in the minimally supersymmetric Lagrangian include a term 

C D -i(V/) T r^° + h.c, 
which is bilinear in the (2-component) fermionic partners 



(Al) 



(A2) 



of the neutral electroweak gauge and Higgs bosons and proportional to the, generally complex and symmetric, neu- 
tralino mass matrix 



/ Mi 
I ° 



-m z sw cp m z sw sp 
y = i •< M 2 mzcwcp -m z c w sp 

I -m z s w cp m z c w cp -\l 

V m z s w sp -mz c w sp -fi 



(A3) 
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Here, Mi, M2, and /1 are the SUSY-breaking bino, wino, and off-diagonal higgsino mass parameters with tan/3 = 
s @l c p = v u/vd being the ratio of the vacuum expectation values v Ut d of the two Higgs doublets, while mz is the SM Z- 
boson mass and sw {cw) is the sine (co-sine) of the electroweak mixing angle 0w After electroweak gauge-symmetry 
breaking and diagonalization of the mass matrix Y, one obtains the neutralino mass eigenstates 

X \ = Niji>l i,i = l,...,4, (A4) 

where N is a unitary matrix satisfying the relation 

N*YN~ 1 — diag (m^o , m^o , m^o , m^o ) . (A5) 

In 4-component notation, the Majorana-fermionic neutralino mass eigenstates can be written as 

*° = (jf )' i = l,... ,A- (A6) 

Their mass eigenvalues m^o can, e.g., be found in analytic form in (56j and can be chosen to be real and non-negative. 
The chargino mass term in the SUSY Lagrangian 

£3 >-)(if)(lr) +li (A7) 

is bilinear in the (2-component) fcrmionic partners 

yjf = {-i^.H^f (A8) 
of the charged electroweak gauge and Higgs bosons and proportional to the, generally complex, chargino mass matrix 

X = ( M % m wV2s P \ (AQ) 

where mw is the mass of the SM W-hosaa.. Its diagonalization leads to the chargino mass eigenstates 

%'-%%' " = 1A (A10) 

where the matrices U and V satisfy the relation 

ITXV- 1 = diag(m-±,m-±). (All) 
In 4-component notation, the Dirac-fcrmionic chargino mass eigenstates can be written as 

*■ = (!)■ (Ai2) 

The mass eigenvalues can be chosen to be real and non-negative and are given by 0] 



m|± 2 = i|M| + 1? + 2m 2 w =p [(M 2 - M 2 ) 2 + 4m^c^ + 4m^(M| + ^ + 2M 2 p s 2/3 )J ' j, (A13) 
while the matrices 

m , TJ - f C+ , if detY >0 ... „ / cos#± sin0± \ 

f7 = O- and V = <^ + , n ' ., , , v ~ n with C± = . c 1 (A14) 

1 03 C + , if det Y < V - sin 6>± cos 6± I x ' 

are determined by the mixing angles 9± with < 9± <-k/2 and 

2v / 2m w (M 2 sin/3 + ^cos/3) 2 V2tow (M 2 cos /3 + fi sin /3) 
tan 20+ = — — -~ t, — o — — and tan 20- = — — -~ - — 5 — — . (A15) 

+ M| - [i 2 + 2m 2 w cos 2/3 M 2 2 - /i 2 - 2™^ cos 2/3 v ; 



1/2 



[1] H. P. Nilles, Phys. Rept. 110 (1984) 1. 



35 



H. E. Haber and G. L. Kane, Phys. Rept. 117 (1985) 75. 

M. Ciuchini, G. Degrassi, P. Gambino and G. F. Giudice, Nucl. Phys. B 534 (1998) 3. 

A. J. Buras, P. Gambino, M. Gorbahn, S. Jager and L. Silvestrini, Phys. Lett. B 500 (2001) 161. 

L. J. Hall and L. Randall, Phys. Rev. Lett. 65 (1990) 2939. 

G. D'Ambrosio, G. F. Giudice, G. Isidori and A. St rumia, Nucl. Phys. B 645 (2002) 155. 
W. Altmannshofer, A. J. Buras and D. Guadagnoli, arXiv:hep-ph/0703200 
N. Cabibbo, Phys. Rev. Lett. 10 (1963) 531. 

M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49 (1973) 652. 
L. J. Hall, V. A. Kostelecky and S. Raby, Nucl. Phys. B 267 (1986) 415. 
J. F. Donoghue, H. P. Nilles and D. Wyler, Phys. Lett. B 128 (1983) 55. 
M. J. Duncan, Nucl. Phys. B 221 (1983) 285. 

A. Bouquet, J. Kaplan and C. A. Savoy, Phys. Lett. B 148 (1984) 69. 
F. Borzumati and A. Masiero, Phys. Rev. Lett. 57 (1986) 961. 
F. Gabbiani and A. Masiero, Nucl. Phys. B 322 (1989) 235. 
P. Brax and C. A. Savoy, Nucl. Phys. B 447 (1995) 227. 
J. S. Hagelin, S. Kelley and T. Tanaka, Nucl. Phys. B 415 (1994) 293. 
F. Gabbiani, E. Gabrielli, A. Masiero and L. Silvestrini, Nucl. Phys. B 477 (1996) 321. 



M. Ciuchini, E. Franco, D. Guadagnoli, V. Lubicz, M. Pierini, V. Porretti and L. Silvestrini, hep-ph/0703204 
W. Beenakker, R. Hopker, M. Spira and P. M. Zerwas, Nucl. Phys. B 492 (1997) 51; 
E. L. Berger, M. Klasen and T. Tait, Phys. Rev. D 59 (1999) 074024. 
H. Baer, B. W. Harris and M. H. Reno, Phys. Rev. D 57 (1998) 5871. 

W. Beenakker, M. Klasen, M. Kramer, T. Plehn, M. Spira and P. M. Zerwas, Phys. Rev. Lett. 83 (1999) 3780. 
E. L. Berger, M. Klasen and T. Tait, Phys. Lett. B 459 (1999) 165; 

E. L. Berger, M. Klasen and T. M. P. Tait, Phys. Rev. D 62 (2000) 095014 [Erratum-ibid. D 67 (2003) 099901]. 
W. Beenakker, M. Kramer, T. Plehn, M. Spira and P. M. Zerwas, Nucl. Phys. B 515 (1998) 3 

E. L. Berger, B. W. Harris, D. E. Kaplan, Z. Sullivan, T. M. P. Tait and C. E. M. Wagner, Phys. Rev. Lett. 86 (2001) 
4231. 

G. Bozzi, B. Fuks and M. Klasen, Phys. Rev. D 72 (2005) 035016. 

S. Dawson, E. Eichten and C. Quigg, Phys. Rev. D 31 (1985) 1581. 

S. Y. Choi, A. Djouadi, H. S. Song and P. M. Zerwas, Eur. Phys. J. C 8 (1999) 669. 

M. Spira, Nucl. Phys. Proc. Suppl. 89 (2000) 222. 

G. J. Gounaris, J. Layssac, P. I. Porfyriadis and F. M. Renard, Phys. Rev. D 70 (2004) 033011. 
A. Bartl, W. Majerotto and W. Porod, Z. Phys. C 64 (1994) 499 [Erratum-ibid. C 68 (1995) 518]. 
M. Obara and N. Oshimo, JHEP 0608 (2006) 054. 
W. M. Yao et al. [Particle Data Group], J. Phys. G 33 (2006) 1. 

M. Ciuchini, A. Masiero, P. Paradisi, L. Silvestrini, S. K. Vempati and O. Vives, hep-ph/0702144 
J. Foster, K. I. Okumura and L. Roszkowski, Phys . Lett. B 641 (20 06) 452. 

T. Hahn, W. Hollik, J. I. Illana and S. Penaranda, |hep-ph/05l"2315"l 

E. Barberio et al. [Heavy Flavor Averaging Group (HFAG)], h ep-ex/0603003| 

A. L. Kagan and M. Neubert, Phys. Rev. D 58 (1998) 094012. 

S. Heinemeyer, W. Hollik, F. Merz and S. Penaranda, Eur. Phys. J. C 37 (2004) 481. 
S. Heinemeyer, D. Stockinger and G. Weiglein, Nucl. Phys. B 690 (2004) 62. 
S. Heinemeyer, D. Stockinger and G. Weiglein, Nucl. Phys. B 699 (2004) 103. 

J. R. Ellis, J. S. Hagelin, D. V. Nanopoulos, K. A. Olive and M. Srednicki, Nucl. Phys. B 238 (1984) 453. 
P. Gondolo, J. Edsjo, P. Ullio, L. Bergstrom, M. Schelke and E. A. Baltz, JCAP 0407 (2004) 008. 
J. Hamann, S. Hannestad, M. S. Sloth and Y. Y. Y. Wong, Phys. Rev. D 75 (2007) 023522. 
J. R. Ellis, K. A. Olive, Y. Santoso and V. C. Spanos, Phys. Lett. B 565 (2003) 176. 
W. Porod, Comput. Phys. Commun. 153 (2003) 275. 

S. Heinemeyer, W. Hollik and G. Weiglein, Comput. Phys. Commun. 124 (2000) 76. 
T. Moroi, Phys. Rev. D 53 (1996) 6565 [Erratum-ibid. D 56 (1997) 4424]. 

B. C. Allanach et al, Eur. Phys. J. C 25 (2002) 113. 

J. A. Aguilar-Saavedra et al., Eur. Phys. J. C 46 (2006) 43. 

M. Battaglia, A. De Roeck, J. R. Ellis, F. Gianotti, K. A. Olive and L. Pape, Eur. Phys. J. C 33 (2004) 273. 
P. Gambino, U. Haisch and M. Misiak, Phys. Rev. Lett. 94 (2005) 061803. 

J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, P. Nadolsky and W. K. Tung, JHEP 0207 (2002) 012. 
M. Beccaria, G. Macorini, L. Panizzi, F. M. Renard and C. Verzegnassi, Phys. Rev. D 74 (2006) 093009. 
B. Fuks, M. Klasen and P. Skands, in preparation. 

M. M. El Kheishen, A. A. Aboshousha and A. A. Shank, Phys. Rev. D 45 (1992) 4345. 



